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Abstract 

The d-dimensional Gaussian free field (GFF), also called the (Eu- 
clidean bosonic) massless free field, is a (i-dimensional-time analog of 
Brownian motion. Just as Brownian motion is the limit of the sim- 
ple random walk (when time and space are appropriately scaled), the 
GFF is the limit of many incrementally varying random functions on 
c?-dimensional grids. We present an overview of the GFF and some of 
the properties that are useful in light of recent connections between 
the GFF and the Schramm-Loewner evolution. 
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1 Introduction 

The d- dimensional Gaussian free field (GFF) is a natural rf-dimensional-time 
analog of Brownian motion. Like Brownian motion, it is a simple random 
object of widespread application and great intrinsic beauty. It plays an im- 
portant role in statistical physics and the theory of random surfaces, partic- 
ularly in the case d — 2. It is also a starting point for many constructions in 
quantum field theory |FFS92l IGK831 IGJ87] . 

The main purpose of this paper is to provide a mostly self-contained math- 
ematical introduction to the GFF for readers familiar with basic probability 
(Gaussian variables, a-algebras, Brownian motion, etc.), but not necessarily 
versed in the language of quantum field theory or conformal field theory. We 
will review the classical continuum constructions (Dirichlet quadratic forms, 
abstract Wiener spaces, Gaussian Hilbert spaces, Schwinger functions, chaos 
decomposition, etc.) and assemble basic facts about discrete Gaussian free 
fields. 

Several results from this paper are cited in a recent work by the author and 
Schramm, which studies contour lines of the discrete Gaussian free field and 
shows that their scaling limits are forms of the Schramm-Loewner evolution 
SLE 4 |SSj . We also expect these facts to be cited in forthcoming work relating 
SLE K to the GFF for other values of k. 

Although [SSj is a long and technical work, it contains an elementary 
twenty-page introduction with many additional references to the history of 
the contour line problem and many other pointers to the physics literature. 
To avoid duplicating this effort, we will not discuss SLE at any length here. 
We also will not discuss the Virasoro algebra or the use of the GFF in the 
Coulomb gas theory (topics discussed at length in several reference texts, 
including |Car90t IDFMS971 IHen99j ). and we generally make no attempt to 
survey the physics literature here. Although this work is primarily a survey, 
we will also present without references several simple results (including the 
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natural coupling of harmonic crystals with the GFF via finite elements and 
the coupling of the GFF and Brownian motion via "field exploration") that 
we have not found articulated in the literature. 

Remark 1.1. In the physics literature, what we call the GFF is often called the 
massless free field or the Euclidean bosonic massless free field — or else intro- 
duced without a title as something like "the field whose action is the Dirichlet 
energy" or "the Gaussian field with point covariances given by Green's func- 
tion." 



2 Gaussian free fields 
2.1 Standard Gaussians 

Consider the space H S (D) of smooth, real- valued functions on M. d that are 
supported on a compact subset of a domain D C M. d (so that, in particular, 
their first derivatives are in L 2 (D)). This space has a Dirichlet inner product 
defined by (/i, = f D (Vfi • V/2) dx. Denote by H(D) the Hilbert space 
completion of H S (D). (The space H(D) is in fact a Sobolev space, some- 
times written Wq(D) or W Q ' 2 (D) |Ada75j .) The quantity (/, /)y is called the 
Dirichlet energy of /. 

Let g be a bijective map from D to another domain D'. If g is a translation 
or an orthogonal rotation, then it is not hard to see that 

/ V(/i o g- 1 ) ■ V(/ 2 o g- 1 ) dx = [ (V/i • V/ 2 ) dx. 

J D' J D 

If g(x) = cx for a constant c, then 

/ V(/i o g- 1 ) ■ V(/ 2 o g- 1 ) dx = c d - 2 [ (V/i ■ V/ 2 ) dx. 

JD' JD 

In the special case d = 2, the equality holds without the c d ~ 2 term. In 
fact, an elementary change of variables calculation implies that this equality 
holds for any conformal map g and any /i,/2 6 H(D). (It is enough to 
verify this for f\, f 2 G H S (D).) In other words, the Dirichlet inner product is 
invariant under conformal transformations when d = 2. (This is one reason 
that the GFF is a useful tool in the study of conformally invariant random 
two dimensional fractals like SLE |SSJ.) 
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When D is a geometric manifold without boundary (e.g., the unit torus 
M. d //L d ), we define H S (D) to be the set of all zero mean smooth functions on 
D, and again we take H(D) to be its completion to a Hilbert space with the 
Dirichlet inner product. 

Note that by integration by parts, (/i,/2)v — (fi,~ A/2)) where A is 
the Dirichlet Laplacian operator and (-, ■) is the standard inner product for 
functions on D. Throughout this paper, we use the notation (fx, / 2 ) v := 
j D (Vfi ■ V/2) dx and (fx, / 2 ) := j D (fxf2) when the integrals clearly make 
sense (even if fx and f'2 do not necessarily belong to H(D) and L 2 (D), re- 
spectively). We also write ||/|| = (f,f)V 2 and ||/|| v = (/, f)f. 

Given any finite-dimensional real vector space V with (positive definite) 
inner product (•,•), denote by /iy the probability measure e~^ v ' v ^ 2 Z~ 1 dv, 
where dv is Lebesgue measure on V and Z is a normalizing constant. The 
following is well known (and easy to prove) | Jan97j : 

Proposition 2.1. Let v be a Lebesgue measurable random variable on V = 
M. d with inner product (•, •) as above. Then the following are equivalent: 

1. v has law /iy. 

2. v has the same law as Ylij=x a j v j where Vx, ■ ■ ■ , Vd are a deterministic 
orthonormal basis for V and the ctj are i.i.d. Gaussian random variables 
with mean zero and variance one. 

3. The characteristic function of v is given by 

Eexp (i(v, t)) = exp(— -||t|| 2 ) 

for all t G M. d . 

4- For each fixed w G V, the inner product (v,w) is a zero mean Gaussian 
random variable with variance (w,w). 

A random variable satisfying one of the equivalent items in Proposition ^. II 
is called a standard Gaussian random variable on V. Roughly speaking, 
the GFF is a standard Gaussian random variable h on H(D). Because H(D) 
is infinite dimensional, some care is required to make this precise. One might 
naively try to define h as a random element of H(D) whose projections 
onto finite dimensional subspaces of H(D) are standard Gaussian random 
variables on those subspaces. However, it is easy see that this is impossible. 
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(Expanded in terms of an orthonormal basis, the individual components of h 
would have to be i.i.d. Gaussians — and hence a.s. the sum of their squares 
would be infinite, implying h G" H(D).) 

We will now review two commonly used (and closely related) ways to 
define standard Gaussian random variables on infinite dimensional Hilbert 
spaces: the abstract Wiener space approach and the Gaussian Hilbert space 
approach. 

2.2 Abstract Wiener spaces 

One way to construct a standard Gaussian random variable h on an infinite 
dimensional Hilbert space H, proposed by Gross in 1967, is to define h as 
a random element not of H but of a larger Banach space B containing H 
as a subspace [Gro67j . To this end, [Gro67j defines a norm | • | on H to be 
measurable if for each e > 0, there is a finite-dimensional subspace E e of H 
for which 

E _L E e /j, e ({x G E : \x\ > e}) < e, 

where is the standard Gaussian measure on E. In particular, we may cite 
the following proposition |Gro67] . (Throughout this subsection, (-, ■) and || ■ || 
denote the inner product and norm of H.) 

Proposition 2.2. If T is a Hilbert Schmidt operator on H (i.e., 

for some orthonormal basis {fj} of H), then the norm \\T ■ || is measurable. 

Write B for the Banach space completion of H under the norm | • |, B' 
for the space of continuous linear functionals on B, and B for the smallest 
a-algebra in which the functionals in B' are measurable. Since each element 
of B' is a continuous linear functional on H, we may view B' as a subset 
of H. Thus B' C H C B. When b G B and / G B', we use the inner 
product notation (/, b) to denote the value of the functional / at b. (When 
/ G H, this is equal to the inner product of / and b in H.) Given any 
finite dimensional subspace E of B' with if-orthonormal basis Vi, . . . , v^, the 
map <pE '■ B — > E given by 4>E{b) = ^2{vj,b)vj is an extension to B of the 
orthogonal projection map from H to E. Let \xe be the standard Gaussian 
measure on E. Gross proved the following: 
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Theorem 2.3. If \ ■ \ is measurable, then there is a unique probability measure 
P on (B, B) for which P((f)^}S) = /j,e(S) for each finite dimensional subspace 
E of B' and each Lebesgue measurable S C E. 

By Proposition 12. 1[ we can restate this as follows: 

Theorem 2.4. If\- \ is measurable, then there is a unique probability measure 
P such that if h is a random variable with probability measure P then for 
any f G B' , the random variable (h, f) is a one-dimensional Gaussian of 
zero mean and variance (f,f) 2 - 

The triple (H, B, P) is called an abstract Wiener space. The example 
that motivated Gross's construction is the standard Wiener space, in which 
H = H ((0, 1)), endowed with the Dirichlet inner product, | • | is the supremum 
norm, and B is the set of continuous functions on [0, 1] that vanish on {0, 1}. 
Using the Hilbert space H(D) (with Dirichlet inner product) we can now 
give a definition: 

Definition 2.5. Given a measurable norm | • | on H(D) and B', B, B as 
above, the Gaussian free field determined by norm | • | is the unique 
5-valued, immeasurable random variable h with the property that for every 
fixed / G B', the random variable (h, /)y is a Gaussian of variance ||/||v- 
Equivalently, h = Y2 a jfj> where ctj are i.i.d. Gaussians of unit variance and 
zero mean and the fj are elements of B which form an orthonormal basis for 
H{D) — and the sum is defined within the space B. (It is not hard to see 
that the partial sums *Y^j=\ a ifi converge almost surely in B [Gro67j .) 

Remark 2.6. We can analogously define the complex Gaussian free field de- 
termined by norm \ ■ | by replacing H, B', and B with their complex analogs 
and writing h = h\ + ih2, where the hi and /12 are independent real Gaussian 
free fields. 

2.3 Choosing a measurable norm 

We now construct one natural family of measurable norms on H(D) using 
the eigenvalues of the Laplacian. Suppose that {ej} are eigenvectors of the 
Dirichlet Laplacian on D which form an orthonormal basis of L 2 (D) endowed 
with the usual inner product and have negative eigenvalues {Xj} (ordered to 
be non-increasing in j). Then an orthonormal basis for H(D) is given by 
fj = (— Aj)~ 1/,2 ej, since integration by parts implies (ej,efc)v = — Ae^) = 
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whenever j ^ k and (fjjjh = {{-\)~ 1/2 Zj, (-Aj) 1/2 ej)v = 1- (This 
choice of the fj is not invariant under conformal transformations of D when 
d=2.) 

The reader may recall that by Weyl's formula, if D C R d is bounded, 
then j 2 / d /(— Xj) tends to a constant as j — » oo. (References and much more 
precise estimates on the growth of Xj are given in |NS05j .) We define powers 
of the negative Dirichlet Laplacian by writing, for each o£l, 

(-ArE^' : =E(-^) a ^' 

a definition which makes sense even when a is not an integer. We then 
formally define C a (D) := (— A) a L 2 (D) to be the set of sums of the form 

for which £ \S *i) "< j e L ( D )- (When a < 0, this sum 
may not converge in L 2 (D), but since (— Xj)~ a is polynomial in j, it always 
converges in the space of distributions on D; see Remark [2 . 8 1 b elow . ) 
Since integration by parts gives 

(f,9h = (f, (-A)y) = ((-A) 1 / 8 /, (-A) 1/2 5f) , 

the map (—A) -1 / 2 gives a Hilbert space isomorphism from L 2 (D) (with the 
L 2 inner product) to H(D) (with the Dirichlet inner product). Thus we may 
write H(D) = £_i /2 (D). 

Similarly, for any a G R, we may view C a (D) as a Hilbert space whose 
inner product (•,•)„ is the pullback of the L 2 inner product, i.e., (f,g) a = 
((-A)-°/,(-A)-°</). We abbreviate ||/|| a := ((-A)- a /, (-A)~ a /) 1/2 for the 
corresponding norm. An equivalent way to define C a {D) is as the Hilbert 
space closure of H S (D) under this norm. 

Proposition 2.7. Suppose D is a bounded domain in R d . Then we have the 
following: 

1. In the space of formal sums ^fijej (or the space of distributions) we 
have C a (D) C Cb(D) whenever a < b. 

2. || • ||b is a measurable norm on C a {D) (where the latter has inner product 
(•, -) a ) whenever a < b — d/4. 

3. When f G C- a (D), the functional g (f,g) is continuous on £ a (D). 
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Proof. The first item is immediate since XX ~ ^j)~ 2a \Pj\ 2 < 00 implies 
^](— Aj) _2b |/9j| 2 < 00. To prove the second, we first write \\f\\b = ||7&_ a /|| a 
where T c := (— A)~ c . Let {gj} be an orthonormal basis for Cb under the 
inner product (•,•)&■ Then ||/||b is a Hilbert Schmidt operator (and hence 
measurable by Proposition 12.21) provided that || T^c/^ || ^ = XX - ^j) 2a ~ 2b < 
00. Weyl's formula implies that this holds provided that ^j 2 ( 2a ~ 2b )/ d < 00) 
which in turn holds whenever 2(2a — 26) jd < — 1, i.e., a < b — d/A. The final 
statement in the proposition is trivial since 

(f,g) = ((-A)*/, (-&) a 9)a = (/, (-A) 2a g) a , 

and (-A) 2a g G C a . I 

Proposition 12.71 implies that, although we cannot construct the GFF h 
as a random element of H(D), we can construct h as a random element of 
B = £ b (D), provided b > ^T 2 -, using the abstract Wiener space definition 
given in Section |2~21 

In particular, when d — 1, we may take b = and define liasa random 
element of L 2 (D). When d — 2, we cannot define ft as a random element 
of L 2 (D) (indeed, from the power series expansion, we expect the L 2 norm 
of h to be almost surely infinite), but we can define h as a random element 
of B = £b{D) for any b > 0. In this case, we may view (h, •) as a random 
continuous linear functional on C-b{D) C L 2 {D) for any b > 0. In general, 
p — > (/i, p) is a random continuous linear functional on C_b{D) whenever 

Remark 2.8. Sometimes it is convenient to restrict attention to smooth, com- 
pactly supported test functions p. Following the usual definition, we say that 
h is a distribution if (h, •) is well defined as a functional on the space H S (D) 
of smooth compactly supported functions and this functional is continuous 
with respect to the topology of uniform convergence of all derivatives. If 
p G H S (D), then (— A) a p G H S (D) C L 2 (D) for each positive integer a, and 
it follows that p G C a (D) for all a. If h G £b(D) for some b, then (h, •) 
is a continuous functional on C C (D) for any negative integer c < —b. This 
implies that the restriction of (h, •) to H S (D) is continuous in the topology 
of uniform convergence of all derivatives (since uniform convergence of all 
derivatives in particular implies convergence in £ C (D)), so h is also a distri- 
bution. Many texts (e.g., |GJ87j ) simply define the GFF to be the random 
distribution determined in this way. Since H S (D) is dense in each of the 
larger spaces Lb(D), we don't lose any information by restricting (h,-) to 
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smooth functions, since there is a unique way to extend (h, •) to a continuous 
function on the larger space. 

Remark 2.9. Let be rotationally symmetric smooth positive bump function 
on M. d whose integral is 1 and which vanishes outside of the unit ball in W d . 
Let f r ,z(x) = r~ d (p((x — z)/r). This function is the density of a probability 
measure on the ball B r (z) of radius r > centered at z G W 1 . Let D r be the 
set of pairs (r, z) for which r > and B r (z) C D. Then the map (r, z) — > / r;2 
is continuous from Z} r to Cb(D) for any & G M. Hence if /i is an instance of the 
GFF defined by one of the norms discussed above, then ip((r,z)) = (h, f rjZ ) 
is a random continuous function from D r to R. Similar arguments show that 
all the derivatives of ip are continuous almost surely. Since the span of the 
f rz is dense, h is almost surely determined by the random smooth function 

2.4 Gaussian Hilbert spaces 

The definition of the GFF in terms of abstract Wiener spaces has an aesthetic 
and practical drawback in that the choice of measurable norm is somewhat 
arbitrary, and it does not yield a description of the random variable (h, f)y 
for general / G H(D). In this section we give a way to make (h, -)v well 
defined as a random variable for each / G H(D) — accepting, of course, the 
fact that / — > (h, /)v cannot be defined as a continuous functional. 

Consider the probability space (Q, J 7 , /x) where Q is the set of real se- 
quences a = {aj},j > 1, J- is the smallest a-algebra in which the coordinate 
projections a — > ccj are measurable, and /i is the probability measure in 
which the oij are i.i.d. Gaussian variables of unit variance and zero mean. 

In the previous section, we defined the Gaussian free field (GFF) to be 
the formal sum h = Yl'jLi a jfj (which converges in a larger space B), where 
the fj are an ordered orthonormal basis for H(D) and the aj are i.i.d. Gaus- 
sians. Now, for any fixed / G H{D) = YlPjfji the inner product (h, /)v is 
a random variable that can be almost surely well defined as the limit of the 
partial sums Y^j=ifij a j- (It i s important here that we fix the order of sum- 
mation in advance, since the sequence fljOij is not necessarily a.s. absolutely 
summable.) 

Now we have a formal definition: 

Definition 2.10. The Gaussian free field derived from the ordered 

orthonormal basis {fj} is the indexed collection Q(D) of random variables 
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(h, /)y described above. 

A more abstract definition, which does not specifically reference a basis 
or an ordering, is as follows. First, we take the following definition from 
|Jan97j : 

Definition 2.11. A Gaussian linear space is a real linear space of random 
variables, defined on an arbitrary probability space (Q,J-',n), such that each 
variable in the space is a centered (i.e., mean zero) Gaussian. A Gaussian 
Hilbert space is a Gaussian linear space which is complete, i.e., a closed 
subspace of L^(Q, J 7 , fi), consisting of centered Gaussian variables, which 
inherits the standard L^(Q, J 7 , /i) inner product: (X, Y) — J XYd/i. We also 
assume that JF is the smallest cr-algebra in which these random variables are 
measurable. 

Note that if Xi, . . . , X n are any real random variables with the property 
that all linear combinations of the Xj are centered Gaussians, then the joint 
law of the Xj is completely determined by the covariances CovXj, XJ = 
E(XjXfc), and it is a linear transformation of the standard normal distribu- 
tion. A similar statement holds for infinite collections of random variables 
[Jan97j . Then we have: 

Definition 2.12. A Gaussian free field is any Gaussian Hilbert space Q(D) 
of random variables denoted by "(h, /)y" — one variable for each / G H(D)— 
that inherits the Dirichlet inner product structure of H(D), i.e., 

E[(M) v (M)v] = M)v- 

In other words, the map from / to the random variable (h, /)y is an inner 
product preserving map from H(D) to G{D). 

By the identity (a, b) = |[(a + b, a + b) — (a, a) — (6, b)], this map is inner 
product preserving if and only if it is norm-preserving — i.e., the variance 
of (h, /)y is (/, /)v f° r eac h / G H(D) — and linear. Thus we have the 
following: 

Proposition 2.13. An H(D) -indexed linear space of random variables de- 
noted (h, f) v is a Gaussian free field if and only if the map from f G H(D) to 
the random variable (h, /)y is linear and each (h, /)y is a centered Gaussian 
with variance (/, /)y. 

Throughout the remainder of this text, we will adopt the Gaussian Hilbert 
space approach and view the variables (h, /)y as being well defined for all 
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/ G H(D). Equivalent ly, we view (h, p) as being well denned for all p G 
(-A)H(D) = C 1/2 . 

Remark 2.14. When p\ and p2 are in H S (D), the covariance of (h,pi) and 
(h, p 2 ) can be written as (— A _1 pi, — A _1 p 2 )v = (— A _1 pi, p 2 ). Since — /Sr l p 
can be written using the Green's function kernel as 

[-A- 1 p](x)= [ G{x,y)p(y)dy, 
Jd 

we may also write: 

Cov[(/i,pi), (h,p 2 )] = / p 1 (x)p 2 {y)G(x,y)dxdy 

JDxD 

When pi = p 2 = p, the above expression has an interpretation in electro- 
statics as the energy of assembly of an electric charge density p (grounded 
at dD), and A _1 p is the electrostatic potential of that density. The energy 
of assembly of a density of charge is the amount of energy required to move 
charge into that configuration starting from a zero-energy configuration (in 
which the potential is everywhere zero). Thus the Laplacian p = (— A)h is, 
at least intuitively, a random electrostatic charge distribution in which the 
probability of p is proportional to 

exp(— energy of assembly of p). 

(See |SSj for more references relevant to this interpretation.) 

2.5 Simple examples 

Let D be the unit torus M. d /Z d . As before H S (D) is the set of smooth func- 
tions on D with zero mean and H(D) is the Hilbert space closure of H S (D) 
using the Dirichlet inner product. An orthonormal basis for the complex 
version of H(D) is given by eigenvectors of the Laplacian, which have the 
form f k (x) = _i e 27r ^- fc , for k G Z d \{0}. Thus, the complex GFF on D 
is a random distribution whose Fourier transform consists of i.i.d. complex 
Gaussians times (27r|A;|) _1 . 

If d > 2, then for any fixed x G D and any fixed ordering of the fc's, the 
partial sums of Y^=i a jfj( x ) diverge almost surely, since the variance of the 
partial sums are given by (27r)~ 2 l^l -2 ; an d this sum diverges when d > 2. 
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When d — 1, the limit h can be defined a.s. at any given x and is a 
complex Gaussian. Since the above sum converges to (27r)~ 2 2£(2) = 1/12, 
the real and imaginary components of h(x) each have variance 1/12. In 
fact, it is not hard to see that the difference h(x) — h(0) can be written as 
(h, f x )sj = (h, 5 X — 5 ) where f x = —A^ 1 (5 X — 5 ) is continuous and linear on 
(0,x) and (a;, 1). By computing dot products of (f x 1 f y )v 1 the reader may 
verify that h has the same law as a multiple of the Brownian bridge on the 
circle, normalized by adding a constant so that it has zero mean. A similar 
argument shows that the one- dimensional GFF on an interval is a multiple of 
the Brownian bridge on that interval, and an even simpler argument shows 
that the one-dimensional GFF on (0, oo) is a Brownian motion. In the latter 
case, we may take [— A~ 1 5 x ](y) = mm{x,y} and G(x,y) = min{x, y}. The 
variance of h(x) = (h,5 x ) is G(x,x) = x. 

Remark 2.15. The definition of the Dirichlet inner product, and hence the 
Gaussian free field, has an obvious analog any manifold on which the Dirich- 
let energy can be defined. In particular, since the Dirichlet inner product 
is conformally invariant when d = 2, the Dirichlet energy has a canonical 
definition for Riemann surfaces. There is also a "free boundary" version in 
which we replace H S (D) by the set of all smooth, mean zero functions on D 
with first derivatives in L 2 (D). 

Remark 2.16. The GFF has a natural dynamic analog in which each (h t , /)v 
is the Ornstein-Uhlenbeck process with zero mean whose stationary distri- 
bution has variance ||/|| v . Thus, instead of taking (h, fj)y to be an i.i.d. 
sequence of random variables, we take (h t , fj)y to be an i.i.d. sequence of 
Ornstein-Uhlenbeck processes parameterized by t. 

2.6 Field averages and the Markov property 

If — Aa = p is constant on an open subset D' C D and equal to zero outside 
of D' (i.e., a is harmonic outside of D'), then we can think of (h, a)y = {h, p) 
as describing (up to a constant multiple) the mean value of h on D' . We will 
retain that interpretation when h is chosen from the Gaussian free field — 
i.e., we think of h as fluctuating so rapidly that it is not necessarily even 
well-defined as a function, but the "average value of h on D'" is well-defined. 

Since Hilbert spaces are self dual, if p is any probability measure on D for 
which / — > pf := J fdp is a continuous linear functional on H(D) (which is 
the case if and only if J2[pfj} 2 < then there is an / for which pg = (/, g)y 
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for all g G H(D), and we have p = -Af G AH(D). 

For example, if d — 2 and p is the uniform measure on a line segment 
L in the interior of D, then ph is well-defined. In this case, the reader may 
check that / — > pf is continuous on H S {D). (The sums X][P./i] 2 < 00 can 
be computed explicitly when D is a rectangle; continuity then follows for 
domains that are subsets of that rectangle.) 

Another important observation is that if Hi and H 2 are any closed or- 
thogonal subspaces of H(D), then (h, -)y restricted to these two subspaces is 
independent. To be precise, denote by Tu 3 the smallest a-algebra in which 
h — > (h, /)v is a measurable function for each / G ii/j. Then it is clear that 
an d ^ff 2 together generate J 7 , and moreover, // is independent on these 
two subalgebras. 

For example, given an open subset U of D, we can write Hu(D) for the 
closure of the set of smooth functions that are supported in a compact subset 
of U. If a G Hjj(D) and b is harmonic in U, then integration by parts implies 
(a, 6)y = (a, —Aft) = 0. Thus Hjj(D) is orthogonal to the closed subspace 
Hjj(D) of functions that are harmonic on [/. 

Theorem 2.17. The spaces H V (D) and H$(D) span H(D). 

Proof. To see this it is enough to show that if / e H S (D), then / can 
be written as a + b, with a G Hjj(D),b G Hjj(D). Roughly speaking, we 
would like to set b to be the unique continuous function which is equal to / 
outside of U and harmonic inside of U, and then write a — f — b. But in 
some cases — e.g., if U is the complement of a discrete set of points — there is 
no b with this property. We will give a slightly modified definition of b and 
show b G H V (D) and / -be Hfr(D). 

Let bg(x) be the expected value of / at the point at which a Brownian 
motion started at x first exits the set Us of points of distance more than 
S from the complement of U. Then as(x) — f — bs(x) is supported on a 
compact subset of U and is clearly in Hu(D). Since Hjj(D) is the closure of 
the union of the Hu s , the a^'s — which are projections onto the increasing (as 
5^0) subspaces H Us — converge to some function a G Hu(D). The b$ thus 
must converge to some b, and b G H^j (since the limit of harmonic functions 
is harmonic) , and f — a + b. I 

For short, we will write T\j = Tu v and = T H ±. The a-algebra is 
one in which random variables of the form (h, /)v = (h, — A/) are measurable 
whenever Af vanishes on U. Intuitively, it allows us to measure the "values" 
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of h outside of U. On the other hand, Tu allows us to measure the "values" 
of h inside of U modulo the harmonic functions on U. The independence of 
the GFF on T\j and can be interpreted as saying that given the values 
of h outside of U, the distribution of the values of h in U is a harmonic 
extension of the values of h on the boundary of U plus an independent GFF 
on U. This property of the GFF is called a Markov property. It holds, in 
particular, if d — 2 and the complement of U is a simple path in D; in this 
case, the .T^-measurable functions measure the values of h along that path 
(or at least the average values of h along subintervals of that path). 

If U is closed and x G D\U, then it is not hard to see that the projec- 
tion f x u of — A -1 ^ onto Hjj has finite Dirichlet energy and that its Lapla- 
cian is supported on the boundary of D\U. Although h(x) = (h, 8 X ) = 
(h, — A~ 1 5 x )y is not a well-defined random variable, we may still intuitively 
interpret (h, f x ,u)v as the "expected" value of h(x) given the values of h 
in U. The reader may check that the function (h, f x ,u)v is almost surely 
harmonic in D\U. (More precisely, since the event u (h, f x ^u)v is harmonic" 
is not in our a-algebra, one shows that the function (h, f x ,u)\? defined on 
dyadic rational points of D\U almost surely extends continuously to a har- 
monic function in all of D\U; this is the same way one proves continuity of 
Brownian motion in, e.g., Chapter 7 of |Dur96j .) We interpret this function 
as the harmonic extension to D\U of the "values" of h on the boundary of 
U. 

2.7 Field exploration: Brownian motion and the GFF 

In this subsection, we describe a simple way of using a space-filling curve 
to give a linear correspondence between the GFF and Brownian motion. 
Roughly speaking, we "explore" the field h along a space-filling curve, and 
the Brownian motion goes up or down depending on whether the values we 
encounter are greater than or less than what we expect. Then each of the 
random variables (h, /)v can be viewed as an appropriate stochastic integral 
of this Brownian motion. Although analogous constructions hold in higher 
dimensions, we will assume for simplicity that d = 2 and D is a simply 
connected bounded domain. 

First, choose / so that A/ is a negative constant on D. Then let 7 : 
[0, 1] — ► D be a continuous space-filling curve. For each t, denote by j t the 
compact set 7QO, £]), and let P t be the projection onto the subspace of 
functions harmonic in D\y t . We also require that 7 remains continuous when 
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it is parameterized in such a way that ||-Pt(/o)||v = * f° r a ^ * e [0> ll/ollvl- 
(This will be the case provided that 7 S is a proper subset of 7 t whenever 
s < 1 In other words, although 7 may intersect itself, it cannot spend an 
entire positive-length interval of time retracing points that have already been 
seen.) 

By decomposing f into its projection onto the complementary subspaces 
^~D\y t an< ^ J~ D \it ) we easily observe the following: 

W(t) :=E((/ i ,/ )v|^\ 7t ) = (h,Pt(fo))v. 

Clearly, W is a martingale, and each W(t) — W(s) is Gaussian with variance 
\s — t\. Hence, W has the same law as a Brownian motion (in the smallest 
a-algebra where each W(t) is measurable). 

Now, the reader may easily verify that the linear span of the functions 
Pt(fo), with t G [0, 1], is dense in H(D). Thus, given the Brownian motion 
W(t), it should be possible, almost surely, to determine (h,fj)y for each j. 
To this end, observe that for any other /, the value 

W f (t) :=E((h,fU^ t ) = (h,P t (f))v 

is also a martingale, and is Brownian motion when time is parameterized by 
Il-f*(/)ll2/v- The question is, how are Wf and W related? 

The answer would be obvious if we had / = aP s (fo) for some fixed con- 
stants a and < s < 1. In this case, Wf(t) = aW(min{s,t}). A similar 
result holds if / is any finite sum of such functions. We have now defined 
Wf(t) for a dense linear space of functions / e H(D), and we may choose an 
orthonormal basis {fj} for H(D) from among that space. Given an arbitrary 
f = Yl a jfj an d any fixed t, we can take Wf(t) to be the limit of the partial 
sums of a jWf. (t) . 

The above discussion gives a linear correspondence between the GFF and 
a Brownian motion. That correspondences of this sort should exist is not 
surprising, given that a Brownian motion can be interpreted as a Gaussian 
free field on the interval / = (0, ||/o|lv) ^ na ^ * s OIir y required to vanish at 
the left endpoint (i.e., H(I) is the Dirichlet inner product Hilbert space 
completion of the set of smooth functions on I that vanish at zero). The 
map that sends the function g s (i) = min{s,t} in H(I) to P s (fo) £ H(D) 
extends to a Hilbert space isomorphism between H(D) and H(I) and the 
correspondence between the GFF and Brownian motion is induced by this 
isomorphism. 
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2.8 Circle averages and thick points 

Fix a domain D C M 2 on which the GFF is defined. For x G D and t 6 R, 
write B x (t) for the mean value of the GFF on the circle of radius e _t centered 
at x (provided t is large enough that the disc enclosed by this circle lies in 
D). The reader may verify that for each x G D, the law of B x (t) is that of 
a multiple of a Brownian motion. If x ^ y, then the GFF Markov property 
implies that the Brownian motions B x (t) and B y (t) grow independently of 
one another once t is large enough so that 2e~* < \x — y\. Let 

A x (t) = ^ „f g(s)e = [°° B x (s)e^ s ds 
J s=t e- S ds J s=t 

be the mean value of the GFF on the disc of radius e~*. 

One can generate various fractal subsets of D by considering the set of 
points x for which B x (t) or A x (t) are in some sense highly atypical. For 
example, one might consider the set of x for which limsup |-B x (t)| < C for a 
constant C, or more generally the set of x for which lim sup | B x (t) — f(t)\ < C 
for a some function /. 

In |HP] the authors fix a constant < a < 2 and define a thick point as 
a point x G D for which ]im A x (t)/t = \fa. They prove that the Hausdorff 
dimension of the set of thick points of a GFF is almost surely equal to 2 — a. 

3 General results for Gaussian Hilbert spaces 

By replacing the Dirichlet inner product with a different bilinear form, it 
is possible to construct different types of Gaussian Hilbert spaces, some of 
which play important roles in constructive quantum field theory [GJ87] . Most 
of the results in the previous section depended heavily on the choice of inner 
product. In this section, we will back up and make some statements about 
the Gaussian free field that are largely independent of this choice. 

3.1 Moments and Schwinger functions 

The moments of the variables (h, p) may be computed explicitly. Suppose 
that pj = — Afj for 1 < j < k. First, we know that E [(/i, p)v] = and 
E [(h, pi)(h, P2)] = (/i,/2)v- We can now cite the following from the first 
chapter of [Jan97j : 
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Theorem 3.1. We have 



fc/2 



E[(h, Pl ) ■ ■ ■ (h, p k )} = HVm* > W v 



M j=l 



where M ranges over the set of all partitions M = {(M^i, Mj^)} of {1, ... ,k} 
into k/2 disjoint pairs. In particular, this value is zero whenever k is odd. 

If each (fj,fk)v is positive (which will be the case, e.g., if each pj is a 
positive probability density function), we can interpret the above value as 
the partition function for a process that chooses a random perfect matching 
M of the complete graph on 1, . . . , k with probability proportional to 



1. e., each edge is weighted by (fi,fj)y. (This is a useful mnemonic, if 
nothing else.) These perfect matchings are simple examples of Feynman 
diagrams (see |Jan97] or [GJ87J). The moment computations above apply 
to general Gaussian Hilbert spaces when (•, -)y is replaced by the appropriate 
covariance inner product. 

Now, for any Xi,...,Xk G D, one would also like to define the point 
moments S/~(xi, . . . ,Xk) := EIX^i h{xj). In the case of the GFF when d > 

2, these h(xj) are not defined as random variables. If they were defined 
as random variables, then, writing p = —A/, we would also expect that 
E(/i, p) = J D Si(x)p(x)dx and more generally 



It turns out that for a broad class of Gaussian and non-Gaussian random 
fields that includes the GFF, there do exist functions (or at least distribu- 
tions) Sk ■ D k — ► M, called Schwinger functions, for which the latter state- 
ment holds, at least when the pj are smooth (Proposition 6.1.4 of [GJ87J). 

In the case of the GFF, the reader may verify that S& is identically zero 
when k is odd and S 2 (xi,X2) = G(xi,x 2 ). For even k > 2, Sk can be 
computed from S 2 using the expansion given in Theorem 13.11 



fc/2 



Y[(-A : pm, .pmJ 
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3.2 Wiener decompositions and Wick products 



From Chapter 2 of jJan97] . we cite the following: 

Theorem 3.2. The set of polynomials in the (h, fj)y is a dense subspace of 
L 2 ({l,F,p). 

The space L 2 (Q, T, /i), endowed with the inner product (X, Y) = K(XY), 
can be viewed as the closure of the direct sum of Hilbert spaces H' n: , each 
of which is the closure of the set of degree n polynomials in (h, fj)x? that are 
orthogonal to all degree n — 1 polynomials (in particular, H :1: = H(D)). The 
decomposition of an element of L 2 (Q, J 7 , p) into these spaces is sometimes 
called the Wiener chaos decomposition (Wiener 1938, Ito 1951, Segal 
1956; see Chapter 2 of |Jan97j ). See Chapters 2 and 3 of [Jan97j for the 
explicit form of the projection operators onto each H' n: , as well as a natural 
orthonormal basis for each H' n: — defined in terms of an orthonormal basis of 
H(D). 

The above implies in particular that if we use the two-dimensional Gaus- 
sian free field to define conformally invariant random sets or random loop en- 
sembles (e.g., via SLE constructions), then any L 2 function of these random 
objects (e.g., the number of loops encircling a given disc) can be expanded 
in terms of this orthonormal basis (although in practice this may be difficult 
to do explicitly). 

Now, given any r)i G H :m: and 772 G H' n: , the Wick product of rji and rj2 
is the projection of 771772 onto H :m+n: . Explicit formulae for the Wick product 
are given in [Jan97j. 

3.3 Other fields 

The "massive" free fields (see Chapter 6 of [GJ87J) may be defined as a 
collection of Gaussian random variables (h, p) with covariances given by 

Cov[(h,pi),(h,p 2 )\ = ((-A + m 2 )- 1 p 1 ,p 2 ), 

where m is a real number called the mass. Here, we can either require the 
test functions p to be smooth (the random distribution interpretation) or 
let them belong to the Hilbert space completion of the smooth functions 
under the inner product ((—A + m 2 ) _1 p 1 ,p 2 ) (the Gaussian Hilbert space 
interpretation). 
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Equivalently, we may consider random variables (h, f)y := (h, /)y + 
m 2 (h, /) where / = (—A + m 2 )~ 1 p. We then have 

Cov[(/i J / 1 )^(/i l / a )^] = (/i,/ a )^ 

The GFF is the case m — 0. Most of the results in this paper have straight- 
forward analogs in the case m^O. 

Among the other fields discussed in quantum field theory (see, e.g., Chap- 
ters 8, 10, and 11 of |GJ87| ) are probability measures that are absolutely con- 
tinuous with respect to the massive or massless free fields and have Radon- 
Nikodym derivatives given by e Z , where V is a "potential" function on 
Q and Z is an appropriate normalizing constant. 

In the more interesting "interacting particle" settings, however, this V is 
undefined or infinite for /x-almost all points in Q. This happens, for example, 
if d > 2 and we write V(h) = j D P(h) where P is an even polynomial. In 
this case, we can define approximations V n to V by writing V n (h) = V(h n ) 
where the h n is a natural approximation to h (e.g., h n could be the partial 
sum Y^j=\ a ifii or ^ could be one of the discrete lattice approximations in 
the subsequent section). We then seek to define a field which is the limit, in 
some sense, of (appropriately normalized versions of) the probability mea- 
sures e~ Vn Z~ 1 fi. When it exists, the limiting measure is in general not abso- 
lutely continuous with respect to /x. See |GJ87j for a mathematically rigorous 
approach to constructing fields, including the so-called P(4>) fields (in par- 
ticular, the celebrated 4 fields), in a way that uses the Gaussian fields as a 
starting point. 

4 Harmonic crystals and discrete approxima- 
tions of the GFF 

4.1 Harmonic crystals and random walks 

Let A be a finite graph with a positive weight function w on its edges. If (pi 
and <p2 are functions on A, we denote their Dirichlet inner product by 

(0i,02)v= w ( e )(Mv) - M x ))(Mv) - M x )) 

e=(x,y) 

where the sum is over all edges e = (x, y) of A. Now, fix a "boundary" <9A, 
which, for now, we can take to be any non-empty subset of the vertices of 
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A. Then the set H(A) of real-valued functions on A whose values are fixed 
to be zero (or some other pre-determined set of boundary values) on dA is a 
|A| — \dA\ dimensional Hilbert space under the Dirichlet inner product. This 
is also a finite dimensional Gaussian Hilbert space, where the probability 
density at is proportional to e - ^^/ 2 . 

When A is a large subset of a rf-dimensional lattice graph L, and dA 
is the set of vertices that border points of L\A, the resulting "discretized 
random surface" model is commonly called the discrete Gaussian free 
field (DGFF), the discrete massless free field, or the harmonic crystal. 

Next, consider a random walk in which each edge e = (x,y) is activated 
by an independent exponential clock with intensity w(e), at which point if 
the position is x, it switches to y, and vice versa. 

Given <fi and x G A\<9A, write Y for the expected value of <f>(x) at the 
first neighbor x of y hit by a random walk starting at y, i.e., 

= Y.e={x,y) W { e ) < t>{ X ) 
£e=(^)™( e ) 

When the weight function w is constant, Y is simply the average value of 
on the neighbors of x; when w is not constant, Y is an appropriate weighted 
average of these values. 

The reader may verify the following one-point discrete Markov property: 
the random variable 4>{y) — Y is independent of the values of <p(x) for x ^ y 
and this random variable has mean zero and variance 1/ Yle={xy) w ( e )- 

This property implies that if we fix the values of on dA (where these 
values are not necessarily equal to zero), then the expected value of any (f>(x), 
for x G A, will be the expected value of 4>(xh), where Xh is the first vertex on 
dA that is hit by a random walk beginning at x. 

We claim that it also implies that Cov[0(x), <f>(y)] is equal to the expected 
amount of time that a particle started at x will spend at y before hitting the 
boundary. (This function in x and y is also called the discrete Green's 
function on A.) To see this, first observe that when y is fixed, both sides 
are discrete harmonic (with respect to the random walk) in A\{y}. It is 
then enough to compute the discrete Laplacian (with respect to the random 
walk) at y itself and observe that it is equal to 1/ J2 e ^ y w ( e ) f° r both sides: 
for the right hand side (the random walk interpretation), this follows from 
well-known properties of exponential clocks. For the left hand side, since Y 
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and 4>(y) — Y are independent and E(0(y) — Y) 2 = 1/ J2 e =(x y ) w i. e )i we have 
E<P(y)(<P(y) - Y) = E(<P(y) - Y) 2 + EY(<P(y) - Y) = 1/ £ w(e). 

e=(x,y) 

As an example, if A is [— C and w — 1, then it is not hard to 
see (from well-known properties of random walks) that the variance of 0(0) 
is asymptotically proportional to n if d — 1, log(n) if d = 2, and a constant 
if <i > 3 (since the random walk is transient in the latter case). 

See jBDGOll IBDZOOl IB1971 IGiaOOj and the references therein for these 
and many more results about harmonic crystals and various generalizations 
of the harmonic crystal. For an analog of the Green's function interpreta- 
tion of variance that applies in the continuum case (known as the Dynkin 
isomorphism theorem), see |Dyn83 , Dyn84a Dyn84b|. 



4.2 Discrete approximations: triangular lattice 

Suppose d = 2, and let L be the standard triangular lattice (the dual of the 
honeycomb lattice). Now, suppose we restrict the GFF on D to the a-algebra 
J r H n where H n (D) is the set of continuous functions that are affine on the 
each of triangles of -L and that vanish on the boundary of D. Since H n (D) 
is a finite dimensional Hilbert space, it is self-dual, and hence a sample from 
the GFF determines an element h n of H n (D), with probability proportional 
to exp[— ||/i n ||y]. Observe also that h n is determined by its values on the 
vertices of the triangular mesh, and that ||/in||v i s equal to ^ \h n (j) — 
h n (i)\ 2 + \ h n {k) — h n {i)\ 2 + \ h n {k) — h n (j)\ 2 where the sum is over all triangles 
(i, j, k) in the mesh. (The area of each triangle is v^3/ (4n 2 ) and the gradient 
squared is §n 2 [|/i n (j) - h n (i)\ 2 + \h n (k) - h n {i)\ 2 + \h n (k) - h n (j)\ 2 }.) Since 
each interior edge of D is contained in two triangles, this is also equal to 
3-1/2 ^ \h n (j) — h n (i)\ 2 , where the sum is taken over nearest neighbor pairs 

('•./)• 

In other words, h n is distributed like 3 1//4 times the harmonic crystal (with 
unit weights) on the set of vertices of -L fl D, where the boundary vertices 
are precisely those that lie on a triangle which is not completely contained 
in D. 

It is not hard to see that the union of the spaces H n (D) is dense in H(D), 
and that for any function / G H(D), f minus its projection onto H n (D) will 
tend to zero in n. Thus, the h n are approximations of h in the sense that 
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for any fixed / G H(D), we have (h n , /)y — * (/i, /)v almost surely. For any 
fixed n, it easy to define contour lines of the continuous function h n . 

4.3 Discrete approximations: other lattices 

Now, again, suppose that d = 2 but that we replace the standard lattice with 
any doubly periodic triangular lattice V (i.e., a doubly periodic planar graph 
in which all faces are triangles). Once again, we can restrict the GFF to the 
functions that are linear on the triangles of \U ■ In this case, the heights 
at the corners turn out to have the same law as a harmonic crystal on this 
triangular lattice graph in which the weight w(e) corresponding to a given 
edge e is given by w(e) = [cot(#i) + cot(#2)]/2, where 8± and 62 are the angles 
opposite e in the two triangles that are incident to e. (Note, of course, that 
the weights are unchanged by constant rescaling.) 

That is (as the reader may verify), the discrete Dirichlet energy 

(K, h n ) v = ^2 w (( x > y))[K{y) - h n (x)} 2 , 

where h n is viewed as a function on the graph, is the same as the continuous 
Dirichlet energy (h n , /i„)v, where h n is extended piecewise linearly to all of 
D. 

In particular, we can get the discrete GFF on a square grid by divid- 
ing each square into a pair of triangles (either direction). In this case, the 
diagonal edges have w(e) = 0, since both 61 and 9 2 are right angles. The 
vertical and horizontal edges have w(e) = 1, since 6\ and 62 are equal to 45 
degrees in this case. (Note that w(e) can be negative if one or both of the 
exceeds 90 degrees. Earlier, we assumed that w was positive; however, since 
the quadratic form corresponding to such a to is still positive definite, the 
definition of the harmonic crystal still makes sense in this setting.) 

Similar results hold when d > 2 if we replace triangles with (/-dimensional 
simplices. 

4.4 Computer simulations of harmonic crystals 

One fortunate feature of the discrete Gaussian free fields is the ease with 
which they can be simulated on computers. Consider first the special case 
that A is an m x n torus grid graph. In this case, the discrete exponential 
functions form an orthogonal basis of the set of complex, mean- zero functions 
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on A, so the elements in a discrete Fourier transform of a discrete GFF on A 
are independently distributed Gaussians with easily computed variances. For 
example, the Mathematica code shown below generates and plots an instance 
of the discrete GFF on an m x n torus. 

ListPlot3D[ Re [Fourier [Table [ 

(InverseErf [2 Random [] -1] +1 InverseErf [2 Random[]-l]) * 

If [j+k==2,0,l/Sqrt[(Sin[(j-l)*Pi/m]~2+ Sin[(k-1)* Pi/n]~2)]], 

{j,m}, {k,n}]]]] 

The code generates a complex GFF and then plots its real part. To parse 
the code, note that the second line simply produces a complex Gaussian 
random variable in Mathematica. For each 1 < j < m, 1 < k < n, the third 
line gives one over the gradient norm of the discrete exponential (x, y) — > 
^(j-i)x^(k-i)y ( w h ere ^ an d £ are mth and nth roots of unity, respectively) 
on the m x n torus (unless j — k — 1, in which case it gives zero). The 
Fourier transform of the corresponding matrix is the Gaussian free field on 
the torus. Now suppose A is a simply connected induced subgraph of the 
torus with values h Q assigned to its boundary; to sample a GFF on A with 
these boundary values, one may first sample an instance h of the GFF on 
the torus and then replace h with h + h where h is the discrete harmonic 
interpolation of the function ho — h (defined on <9A) to all of A. The most 
time-consuming part of this algorithm is computing the discrete harmonic 
interpolation (but it is not hard to compute an approximate interpolation). 

4.5 Central limit theorems for random surfaces 

Kenyon in [KenOlj proved that random domino tiling height functions (in 
regions with certain kinds of boundary conditions) converge to the GFF as 
the mesh size gets finer. (To be precise, |Ken01j studies random discrete 
height functions h e — for which the lattice spacing is e — and shows that for 
smooth density functions p, (h t , p) converges in law to (h, p) = {h, — A~ 1 p)y.) 
Also, |NS97| IGOS01] give a similar Gaussian free field convergence result for 
a class of discretized random surfaces known as Ginzburg-Landau V0 
random surfaces or anharmonic crystals and [GOSOlJ shows further 
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20 

Figure 4.1: Discrete Gaussian free field on 20 by 20 grid with zero boundary 
conditions. 

that certain time- varying versions of these processes converge to the dynamic 
GFF. 

See |She05] for more references on random surface models (assuming both 
discrete and continuous height values) with convex nearest neighbor potential 
functions and discussion of the division of the so-called gradient Gibbs 
measures into smooth phases (in which the variance of the height difference 
between two points is bounded, independently of the distance between those 
points) and rough phases (in which the variance of these height differences 
tends to infinity as the points get further apart). An important open question 
is whether every two-dimensional rough phase has a scaling limit given by a 
linear transformation of the GFF. 
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